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Abstract 

We propose a new interpolating field for Sii(1535) to determine its mass 
from QCD sum rules. In the nonrelativistic limit, this interpolating field 
dominantly reduces to two quarks in the s-wave state and one quark in the p- 
wave state. An optimization procedure, which makes use of a duality relation, 
yields the interpolating field which overlaps strongly with the negative-parity 
baryon and at the same time does not couple at all to the low lying positive- 
parity baryon. Using this interpolating field and applying the conventional 
QCD sum rule analysis, we find that the mass of Su is reasonably close to the 
experimentally known value, even though the precise determination depends 
on the poorly known quark-gluon condensate. Hence our interpolating field 
can be used to investigate the spectral properties of Sii(1535). 
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The negative parity nucleon resonance Sii(1535) has always been an object of study 
when constructing models of the baryon. A good fit to its mass and decay widths provides 
further test of the validity of the model, which in return provides a deeper understanding of 
the structure of the resonance itself. 

Recently, in relation to the present and planned experiments of r^-production on a nu- 
cleon from electromagnetic or hadronic probes at MAMI (Mainz) [Q, ELSA (Bonn) [0] and 
CEBAF, there has been a renewed interest in Sn. Because the mass of Sn is just above the 
r/N threshold, 77-pro duct ions in these experiments are dominated by s-channel Sn resonance 
contribution. As a consequence, a detailed study of its properties, such as the 77-nucleon-Sii 
coupling and transition form factors, can provide important constraints in the analysis of 
these experimental data. So far, most theoretical works P] have been focused on repro- 
ducing the experimental data by using the empirical parameters related to the Sn, which 
have been explicitly introduced in their effective models. To make a connection between 
these effective models with presumably the fundamental theory of QCD, it is important to 
understand how the spectral properties of Sn are generated from QCD order parameters. In 
this way, phenomenological aspects of Sn can be predicted rather model-independently. In 
this letter, as a starting point for this kind of analysis, we will introduce a new interpolating 
field for the Sn suitable for conventional QCD sum rule analysis. 

The QCD sum rule @] is a useful tool to study the spectral properties of the hadrons. It 
has been widely used in a number of occasions; in calculating properties of hadrons, coupling 
strength between hadrons, form factors and recently medium dependent properties . 
The starting point of the QCD sum rule is to introduce an appropriate interpolating field 
for the hadron of concern. In general, the interpolating field is constructed from quark and 
gluon fields, whose specific form is usually determined by considering the hadron's quan- 
tum numbers and its nonrelativistic quark wave functions. Then the two-point correlation 
function of such interpolating field is introduced and calculated in the operator product 
expansion. By matching the resulting correlation function with its phenomenological coun- 
terpart, one can make prediction on the hadron's spectral properties in terms of the QCD 
order parameters. Clearly, the success of the prediction is highly dependent upon the choice 
of the interpolating field. 

Recently, Jido Kodama and Oka P] have proposed a technique to separate out the 
contribution of the negative-parity baryon from QCD sum rules using the usual interpolating 
fields of the nucleon. They succeeded in obtaining a reasonable value for the mass of Sn by 
adjusting the QCD input parameters and at the same time optimizing a linear combination 
of the two independent nucleon interpolating fields. However, there one has to take the 
difference between the two independent sum rules of the nucleon each proportional to the 
Dirac structure of 1 and ^ so as to subtract out the nucleon contribution, which has a large 
overlap with the current. This fact makes it difficult to investigate other properties, such 
as the form factors or the couplings, of the Sn resonance within the conventional sum rule 
approach. The problem is that the two independent nucleon currents used so far either 
couple strongly or weakly to both the nucleon and the low-lying resonances at the same 
time, such that it is difficult to select out only the negative parity part of the resonances in 
either case. To overcome this difficulty, we will construct a current that couples weakly to 
the nucleon and strongly to the Sn. It should be noted that the reason why the nucleon 
couples strongly with the commonly used loffe current is that, this current has a strong 
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overlap with the nonrelativistic quark wave function of the nucleons in which all quarks are 
in the s-wave state. Therefore, we will construct a current whose nonrelativistic limit has 
a strong overlap with two quarks in the s-wave state and one quark in the p-wave state. 
This is the picture for Sn both in the bag model or nonrelativistic quark model. This 
implies that we have to go beyond the two independent nucleon sum rule and introduce an 
appropriate covariant derivative in the current, for which there are many choices. We will 
choose a minimal approach and construct our current so that it does not couple to the low 
lying positive-parity nucleon but couples strongly to the negative-parity baryon. As will 
be shown by QCD sum rule analysis, our current gives a reasonably stable plateau for its 
mass in the Borel curve. Using this current, it is possible to apply conventional sum rule 
approaches to investigate other properties of the Sn resonance, which will be a subject of 
future study 

First, we start with a short discussion of the proton interpolating field. The interpolating 
field of the proton is constructed from two u-quarks and one d-quark by assuming that all 
three quarks are in the s-wave state. Specifically, one up and the down quarks are combined 
into an isoscalar diquark and the other up quark is attached to the diquark so that the 
quantum numbers of the proton are carried by the attached up quark. The proton should 
be a color singlet and its interpolating field should transform as the spinor under the parity 
transformation. Then, it is found that there are two interpolating fields possible ITOHl^, 



Vl = ^abc{ulCdb)'^^Uc , (1) 

ri2 = eabc{uJC'y5db)uc , (2) 

where r]2 contains the valence quark wave function in the nonrelativistic limit. In general, 
the interpolating field for a proton is an arbitrary linear combination of these fields, 

ri{t) = 2eabc[ t (u^Cdbh^Uc + (uJC754)mc] • (3) 

For t equal to —1, the resulting interpolating field can be shown to have no direct instan- 
ton contribution and, at the same time couple strongly to the chiral symmetry breaking 
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Sii is the lowest resonance of the nucleon with negative parity. As suggested by the 
bag model |T^ or nonrelativistic quark model |]T^, one quark in Sn is believed to be in 



the p-wave state relative to the other two quarks. We want to construct a current, which 
in the nonrelativistic limit overlaps largely with this quark field configuration. One way to 
accomplish this is to put a covariant derivative z ■ D to one of the quarks in Eq. (j^). Here, 
we have introduced a four vector 2;^ such that z ■ q = and z"^ = —1. The property of z^ is 
chosen to make the covariant derivative orthogonal to the four vector q carried by the Sn 
so that in the rest frame, z ■ D reduces to the derivative in the space direction. Of course, 
there are three possible choices for putting in the derivative inside the nucleon current. In 
this work, we take the interpolating field with the covariant derivative acting on the d-quark 
instead of other two u-quarks : 

r]N- (t) = 2eabc[ t iu^C{z ■ D)db)i^u^ + {ulC-i^{z ■ D)db)uc] . (4) 

This choice is preferred because the resulting correlator gets nonzero contribution from the 
lowest-order chiral breaking term, (gg), which is important in the nucleon sum rule. 
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With this current, we now move on to calculate the time-ordered correlation function 
defined as 



U{q) = J dxV-^z (0|T[r/^-(x)?7^-(0)]|0) , (5) 

where |0) denotes the QCD vacuum. Implementing the properties of 2^, it is easy to show 
that the correlation function should have the following form, 

n{q) = n,{q',z') + n,{q^z') 4. (6) 

Calculation of these two scalar functions is carried out readily in the operator product 
expansion. After some manipulations, we obtain up to dimension 8 operators. 
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{m){9sqo-QqW + t) , (7) 



nr(g) = ^^^(^^ - mm) - ^^^(3^^ + 2t - ^s-qo ■ Qq) 

(8) 

Note that, in calculating these, we have made use of the properties, z ■ q = ^ and z^ = — 1. 
Few remarks are in order. The Wilson coefficients of the dimension six condensate {{qq)"^) 
in n°P'^(g) can be shown to be zero. We added the contribution from the dimension eight 
operator {{qq){gsqc^ ' Gq))^ they contribute in tree graph and hence expected to be the 
important power correction for n°P*'. For ni(g^,2;^), the leading contribution comes from 
dimension three condensate, {qq)., which constitutes the important chiral symmetry breaking 
operator in predicting the nucleon mass in the nucleon sum rule. In addition, another chiral 
breaking operator, {gsqty ■ Qq), contributes to the Sn sum rule. This operator does not enter 
to the usual nucleon sum rule where loffe current is employed. However, as will be discussed 
below, our prediction for Sn mass is sensitive to the value of this dimension-five condensate. 
Note also that the factorization hypothesis has been employed in the calculation of higher 
dimensional operators. We did not include the three-gluon condensate {Q^) which is an 
order of a^J"^. 

After Borel transformation, Eqs. (0)and (|^) become, 

WJ>\M) = -(21^2 + m + 21) + -(45^2 + lot + A5)(—g^) 

1 ^ ' 29x3x5x7r4^ ' 2^ x 9 x ^ 'n' 

+^{t^ + t){qq){9sq(7-gq) , (9) 

nr(M) = -i^,{qq)if - 1) + 1^2^9sqo ■ Qq) (3t^ + 2t - 5) 

'^g-){qq){e-l), (10) 



192 ' TT 

where M denotes the Borel mass 
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For the phenomenological side of the sum rules, first we consider the matrix element of the 
interpolating field [Eq. (Q)] between the QCD vacuum and a baryon state with momentum 
q and spin s. Due to the specific form of the field, we should have 

{0\vmQ,s)=^XBu{q,s) , (11) 

where u{q, s) is a Dirac spinor and is the strength with which the baryon couples to the 
interpolating field. Normally, the interpolating field couples not only to the positive-parity 
baryon but also to the negative-parity baryon. This implies that, when combined with 
Eq. (pJ]), the phenomenological form of the correlator should be 

Here we retain the lowest resonance of positive parity which we identify as the nucleon(A^), 
and the lowest resonance of negative parity which we identify as Sn {N~). We have put 
all other resonances with higher masses into two separate continua, one for the positive- 
parity resonances which we denote by s+, the other for the negative-parity resonances which 
we denote by s_. Xn and X^- represent the strengths with which our interpolating field 
couples to nucleon and Sn, respectively while the corresponding masses are denoted by Mjsf 
and M]y-. 

In general, the couplings, Xn and A^r- are functions of t. To obtain the sum rule for 
Sii, the parameter t need to be chosen such a way that the resulting interpolating field 
does not couple to the positive-parity baryons while strongly couples to Sn. This can be 
achieved by applying the finite energy sum rule as described in Ref. We will come back 
to this procedure later. However, once t is chosen this way, there is no contribution from the 
nucleon, and therefore one can derive the sum rule for Sn in the conventional way. Following 
the usual steps to include the continuum [see for example Ref. |]I2|] and taking the Borel 



transformation of the resulting expressions, we get the phenomenological side of the sum 
rule assuming two different continuum thresholds, 

nf-(M) = X%.e-K~/''^' + (21t!+10t + 21) 1 s « 
1 ^ ^ ^ 29 X 3 X 5 X 7r4 2^ 

Hf ^'^(M) = Xl-M^-e-""'^-/'''' - M(t^ - l)^(Mf + M^) 
_ (^,|_^^3^, + 2t - 5)i(Ml + Ml) 

+ Ml) , (13) 



where we have defined. 

Ml = M^e-'^"'" , (14) 
s± ■ 
M2- 



Ml = M\l + ^)e-^±/*^^ , (15) 
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M, 



M. 



M» 1 + ^ + 



+ 



„3 



-s±/M^ 



M2 2M4 6M6 



(16) 
(17) 



By equating these phenomenological side with the OPE side of the sum rule [Eqs. 
(p!0|)] and taking the ratio, we obtain the sum rule for M^v-, 



7V- 



X 



A(Mi)+A(M6 



if-l){qq) + 



A{Mt: 
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A(Mf ) , 

—-^(3f + 2t-5){gsqa-gq) 



A{Ml)+A{M^ 
24 

A{Ml) + A{Mt 



9607r2 
A{Mt)+A{Mt 



256 X 9 



-(2ir + 10t + 21) 
-(— 6^') (45^2 + lot + 45) 



TT 



47r^ 



{t' + t){qq){9sqa-gq) 



-1 



(18) 



where we have defined 

(19) 
(20) 
(21) 
(22) 

Now, the QCD parameters appearing in Eq. (|18]) should be properly chosen to predict 
Mjv-. The quark condensate can be deduced from Gell-Mann— Oakes— Renner relation 
combined with the current quark masses estimated from chiral perturbation theory fT^ : 





= M^ 


-Ml 




= M^ 


-Ml 


A{Ml) 


= 


-Ml 


A{Ml) 


= M' 


-Ml 



{qq) = -(230 ±20 MeV)^ . 



(23) 



The gluon condensate can be estimated either from leptonic decays of vector mesons or from 
the charmonium spectrum [Q: 



{^g^) = (350 ±20 MeV)^ . 

TT 



(24) 



The most important parameter for our prediction is the quark-gluon condensate which 
is usually expressed in terms of the quark condensate {qq) : 



{9sq(y-gq) =2Xl{qq) 



(25) 



Here, the parameter A^, which represents the average virtual momentum of vacuum quarks, 
is not precisely known. It has been estimated from a number of studies; = 0.4 ±0.1 GeV^ 
from the standard QCD sum rule estimate |T^, A^ = 0.7 GeV^ from QCD sum-rule analysis 
of the pion form factor [16|, A^ = 0.55 GeV^ from lattice calculations fl^, and somewhat 
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larger value from the instanton liquid model. In our discussion, we will investigate the 
sensitivity of our prediction by changing between 0.4 — 1 GeV^. 

Given the QCD parameters, we need to choose an optimal t such that the interpolating 
field does not couple to the positive-parity nucleon at all while maximally couples to the 
negative-parity nucleon. To do this, we use the sum rule derived by JKO in Ref. where 
they have used the "old-fashioned" correlation function defined by 

n(g) = 1 1 d''xe'i^e{x,){Q\m-{^)vN-m^) ■ (26) 



This function is analytic in the upper-half region of the complex go plane. Therefore, using 
the nucleon current and assuming, 

Im n(go, q = 0) = 7o^(go) + B{q,), (27) 
one can derive the following sum rules for a large energy Q, 

r dqo[B^P%qo) - B^"^^^ {qo)]W (qo) = , (28) 
Jo 

where W{qo) is the weighting function. The advantage of using the "old-fashioned" sum 
rule is that it allows us to determine the optimal t independent of baryon masses which we 
want to calculate ultimately. 

As stated above [Eq.(0)], we assume that our phenomenological side has explicit con- 
tributions from the nucleon (A^) and the Sii{N~). That is, we have 

-Im [nP^^'^'^(go)] = \%^^^Siqo - M^) + A^„^^<5(go - M^-) + continuum . (29) 



Then we derive the following sum rule from Eq.(pq) with Vr(go) = 1 

A(s+) - B{s^) = Xl , (30) 
A{s.) + B{s.) = Xl- , (31) 

where in our case, 

A{s^)= [^dqo\-^—^ j(21t2 + 10t + 21) 

+ ^nx9x7r ^^^^' + lOt + A5){^g') + -n6{qo){qq){gsqa ■ gq){t' + t)] , (32) 
B{s^) = [^dqo\- ^{f - l){qq) + + 2t - 5)(^7,ga ■ Gq) 

Jo L Z^TT Z'TC 

■K{t^ - 1) ^Q!„ 2\/- 



192 ^^^'m- (33) 

Here we fix the positive-parity continuum threshold to the second lowest resonance of its 
kind, s+ = (1.44 GeV)^, and for the negative-parity threshold, we take s_ = (1.65 GeV)^ as 
it is the next higher resonance with the negative-parity. 
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Eq.(|30|) is a kind of duality relation between the quark and nucleon state. It equates 
the integrated positive-parity spectral strength within the duality interval (0 ~ y/s+) of 
the OPE part to that of the hadronic counterpart, which is the coupling strength of the 
nucleon. Note that the right hand sides of Eqs. (|30|) , (|3TD are positive definite. Now, we 
look for the optimal t which makes the left hand side of Eq. ( pO]) zero while maximizing A^- 
from Eq. (pi]). Our expectation is that such t yields the interpolating field which couples 
strongly to the negative-parity baryon, specifically Sn. Since the quark-gluon condensate is 
not well known, we determine t for the following quark-gluon parameter, 

= 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1 GeV^ . 

These values are arbitrary chosen from the acceptable range of A^ between 0.4 — 1 GeV^. 
[See Eq. ( ^Sf ) for the definition of A^.] With the following parameter set. 



(qq) = -(0.23 GeV)^ (34) 
{^g^) = (0.35 GeV)^ , (35) 



TT 



we obtain the optimal t for each A^ which are listed in Table |. Note that our result is not 
sensitive to the quark or gluon condensate within the error bars as indicated in Eqs. (^31) 
and (p4D. In general, there are two values of t which make the left-hand side of Eq. (|30D 
zero. We choose t which yields larger value of A^_ in Eq. (31). 



With this optimal t value, we are now in a position to calculate Mj^-. Substituting 
all the parameters into Eq. (^) we obtain the Borel curve for Mj^-. Some of them are 
shown in Fig. |1]. The solid line is for A^ = 0.4 GeV^, the dashed line is for A^ = 0.6 GeV^, 
and the dot-dashed line is for A^ = 1 GeV^ . Note that all three curves have the stable 
region from which the spectral mass can be obtained. At A^ = 0.6 GeV^, our prediction for 
M^v- obtained from the stable Borel plateau is 1.59 GeV. Similarly, at A^ = 1 GeV^, our 
prediction is M^r- = 1.42 GeV. All other predictions for M^r- with different A^ values lie 
within these two limits as shown in Table |I[ The best prediction for Sn mass is obtained 
for Ag = 0.435 GeV^ or A^ = 0.69 GeV^. However, for A^ = 1 GeV^, our prediction is off by 
7% from the experimentally known Sn mass, 1.535 GeV. To make more reliable prediction, 
it is important to narrow down the range of the quark-gluon condensate value from other 
studies. However, it should be emphasized that our prediction is obscured only by 7% from 
the barely known quark-gluon condensate. That is, for a given value of A^, the optimal 
procedure that we have employed chooses the interpolating field which couples very strongly 
to Sii. 

As we have succeeded in predicting the Sn at least in qualitative level, it would be 
interesting to understand the origin of the large mass splitting between Sn and nucleon. 
It should be noted that the dominant portion of our prediction comes from the dimension- 
five quark-gluon operator. The contribution from this dimension-five operator is enhanced 
by the covariant derivative acted on the d-quark of the nucleon interpolating field because 
the covariant derivative modifies the quark configuration to increase the average virtual 
momentum of the current quark. However, this is not the case for the nucleon sum rule 
where the quark configurations are mostly in the s-wave state. Indeed, with lojfe^s choice for 
the nucleon interpolating field, the quark-gluon condensate does not enter in the prediction 
of nucleon mass. Therefore, within our picture, the large mass splitting between nucleon 
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and Sii is mainly driven by the enhancement of the quark-gluon condensate which increases 
the contributions from the chiral-symmetry breaking terms. 

Another interesting feature of our result is that, for our current, we can always find 
suitable linear combination of currents such that it has a strong coupling with the negative- 
parity baryon and zero coupling to the positive-parity baryon. This justifies the use of 
Eq.(|l8]) or any other sum rule constructed with this current under the assumption that the 
lowest mass resonance that couple to this current is the Sn. This is not always possible 
for any other interpolating field for the nucleon, especially when the current contains no 
derivative. For example, using Eq.(|^) as an interpolating field, it is not possible to find 
any value of t such that the current has a strong coupling to the negative-parity baryon 
and zero coupling to the positive-parity baryon. Therefore, our interpolating field can be 
used to construct any conventional sum rule to study other properties of the Sn relevant in 
reactions involving electromagnetic or hadronic probes 

In summary, we have introduced a new interpolating field for Sn motivated by the bag 
model or nonrelativistic quark model. Using the interpolating field, we have calculated the 
correlation function in the operator product expansion up to dimension eight. To predict 
Sii mass, we have optimized the linear combination at a given QCD parameter while 
fixing other well known parameters. A good prediction for Sn mass is obtained under this 
formalism and its value contains at most 7% error due to the poorly known QCD parameter 
Xg. The best value for Sn mass is obtained for = 0.435 or 0.69 GeV^. 
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TABLES 



TABLE I. Our prediction for Sn at given with optimally chosen t. The listed Mjv- is 
obtained from the stable Borel plateau. 



(GeV^) Optimal t M^- (GeV) 

04 497 L49 

0.5 -47.52 1.59 

0.6 -4.59 1.59 

0.7 -2.91 1.53 

0.8 -2.40 1.48 

0.9 -2.21 1.44 

1.0 -2.04 1.42 
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FIGURES 

FIG. 1. Our prediction for Sn mass versus Borel mass. The solid curve is for = 0.4 GeV^, 
t = 4.97 and the dashed curve is for Xl = 0.6 GeV^t = -4.59. The dot-dashed curve is for 
A2 = 1 GeV^ t = -2.04. 
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